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Pair production of light top squarks at the 8-TeV LHC can be used to probe the gaugino-Higgsino 
sector of the Minimal Supersymmetric Standard Model. The case where the lightest neutralino is 
a mixture of Bino and Higgsino, satisfying the thermal dark matter relic density, is investigated. 
In such a scenario, the lightest top squark decays mostly into (i) a top quark plus the second or 
third lightest neutralino, and {ii) a bottom quark plus the lightest chargino, instead of a decay 
scenario of the lightest top squark into a top quark and the lightest neutralino. Final states with 
> 2 jets, dileptons, and missing energy are expected in a subsequent decay of the second or third 
lightest neutralinos into the lightest neutralino via an intermediate slepton ("light sleptons" case) 
or Z boson ("heavy sleptons" case). The opposite-sign same flavor dilepton mass distribution after 
subtracting the opposite-sign different flavor distribution shows a clear edge in the case of light 
sleptons. The significance for discovering such a scenario is calculated with optimized cuts in both 
light and heavy sleptons cases. 



I. INTRODUCTION 



Results from the 8-TeV LHC (LHC8) have put bounds 
on the masses of the gluino [g) and first two generation 
squarks (g) in models of supersymmetry (SUSY). If they 
have comparable masses, the exclusion limits reach to 
approximately 1.5 TeV at 95% CL with 13 fb"! of in- 
tegrated luminosity [l]-[3]. The third-generation squarks 
have also actively been searched, putting experimental 
bounds on parameter space in certain decay modes [HIS]. 

Light top squark searches so far have been carried out 
on the case where the lightest neutralino {x\) is mainly 
a Bino and the second lightest neutralino (^2) mainly a 
Wino. In such a scenario, the lightest top squark (ti) 
decays to and a top {t) quark at a branching fraction 
{B) of nearly 100%. In a recent paper by some of the 
current authors [7], the trijet invariant mass MS was used 
to reconstruct top quarks in fully hadronic final state of 
events with at least four non-& jets, at least two h jets and 
large missing energy. There have been other approaches 
to probe this decay by different groups [8j. Reference [5] 
has studied the it decay in the scenario where Xi and the 
lightest chargino (x^) are purely Higgsino. 

Due to small electroweak (EW) production the bounds 
on the neutralinos and charginos are much weaker. This 
sector, along with the sleptons, plays a crucial role in the 
dark matter physics of supersymmetric models. In the 
i?-parity conserving Minimal Supersymmetric Standard 
Model (MSSM), x? is typically the dark matter (DM) 
candidate. If Xi is purely a Bino, its relic density tends 
to be large since the annihilation cross-section is smaller 
than the required thermal annhilation rate 3 x 10~^^ 



OTP? /sec} . One way to obtain the correct relic density is 
to consider a thermal, well-tempered xj which is a mix- 
ture of Bino and Higgsino O , while having X2 and 
X3 as primarily Higgsinos.^ 

The purpose of this paper is to utilize ti decay to probe 
the dark matter sector in a scenario with Xi as a Bino- 
Higgsino mixture and X2 3 as mainly Higgsinos. All three 
are lighter than the lightest top squark, which is in the 
sub- TeV range. The main theoretical motivation for con- 
sidering a light top squark as well as light Higgsinos is 
Naturalness, while the motivation for the presence of a 
light Bino is to obtain the correct relic density for Xij 
since if a sub- TeV Xi is purely Higgsino, the relic density 
is too small P7] . 

In such a scenario, the lightest top squark mainly de- 
cays into tx2,3 and 6x? , followed by X2,3 ^ ^Xi or ^'Xi 
(via an intermdiate slepton state) and x^ ~^ ''^Xi- The 
final state in tit\ events has dileptons with jets and miss- 
ing energy {pr)- The cases of light slepton (whose mass 
is between x? and X2) and heavy sleptons are considered 
in events with at least two leptons, jets and fjj:. 

The dilepton final states investigated in this paper can 
lead to a quite robust ii search. The cross section for 
1x1*1 production is appreciable at the LHC8 for the mass 



^ The correct relic density for pure Bino can still be obtained by 
different methods, such as coannihilation 1101 . slepton exchange 
for lighter Bino masses, resonance, or non-thermal production 
without further annihilation 1111 . 

^ Other options include non-thermal Winos or Higgsinos |14) . 
multi-component dark matter |15l . or well-tempering with Bino 
and Wino. We note that if is purely a Wino, the annihilation 
into W+W- final states is in tension with Fermi data for Wino 
mass below - 250 GeV [16]. 
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range between 300 and 700 GeV. It is shown that the 
SUSY combinatoric and SM backgrounds are reduced 
by performing a opposite-sign same flavor (OSSF) mi- 
nus opposite-sign different flavor (OSDF) subtraction. If 
slepton masses are between X2 Xi; an edge in the 
dilcpton mass distribution could be visible due to higher 
branching fractions of X2 3 ~^ '^Xi decays. The presence 
of a fe-tagged jet in the final state is key to inferring the 
production of a third-generation squark. 

The most dominant background for titl events is ti 
plus n-jets processes. ATLAS has reported bounds on 
the m^-rrij^o plane for single lepton, dilepton, and fully 
hadronic final states using 13 fb^^ of data at y/s = 8 TeV 
[HIS]. Bounds have also been reported for the ii — > bxf 
decay mode. 

It is worthwhile to point out that for probing the dark 
matter sector with mostly-Bino Xi and mostly- Wino 
vector boson fusion processes are a powerful tool [IB]. 
This is somewhat complementary to the strategy in this 
paper. 

The rest of the paper is structured as follows. In Sec- 
tion |llj the results for the "light slepton case" m^^ > 

m^o ,~± > nij > m^a are described, while results for the 

X2JA1 * ^1 

"heavy slepton case" m? > m*, > m~o ,-± > are in 

1 X2'Al Xl 

Section [ni] The dark matter relic density for the differ- 
ent scenarios are discussed in Section [TV] followed by our 
conclusions in Section IVl 

II. SCENARIOS WITH LIGHT SLEPTONS 

In this section, scenarios which satisfy the following 
mass relation are studied: 

mf > m^o ,m^o ,m~± > mr > m^o . (1) 

The possible ii decay modes are: 

ii ^ t x! (2) 

ii txl^t /=F[(*)± ^ t l^l^xl, (3) 

h ^ bxt^b Wxl (or b qq'xl) (4) 

ii ^ bxt^b Zx+ (5) 

Throughout this paper inclusion of charge conjugate 
modes is implied. The last mode is allowed when the 
X2 is lighter than ii. 

It is clear that one obtains an edge in the dilepton 
invariant mass distribution as well as Z-peak depending 
on the size of B{ii bxf) value. A mass spectrum at 
our benchmark point is displayed in Table |l] 

The mass spectrum of the model is determined using 
ISASUGRA [2^. The spectrum is then fed to PYTHIA [21], 
which generates the Monte Carlo hard scattering events 
and hadron cascade. These events are passed to the de- 
tector simulator PGS4 [52]. We use darkSUSY [23] for rehc 
density computations. 

The latest ATLAS exclusion limits on the chargino- 
neutralino plane are given in |24j , in events with 31 + f/j 



in 13.0 fb^^ of data at y/s — 8 TeV. Exclusion plots 
for the direct production of charginos, neutralinos, and 
sleptons are also given in ^25; , in final states with at least 
two hadronically decaying rs and pr- Our benchmark 
values of Xi, xj are outside the exclusion limits given by 
ATLAS. 

TABLE I: SUSY masses (in GeV) at "light slepton" bench- 
mark point. 



Particle Mass (GeV) B 



ii 


500 


17% (tx§), 22% {txD, 8% (tx?) 
53% {hxt) 




175 


100% {fl) 




176 


88% {li) 


X? 


164 


22% {IvxD 


I 


144 


100% Hxl) 


X? 


113 





In the benchmark scenario, the mass difference be- 
tween X2.3 and Xi is around 63 GeV, and thus an edge 
in the dilepton invariant mass distribution may be ex- 
pected around this value. The final state of 2 jets + 2 
leptons -|- fjT events arises mostly from a combination of 

the <i bxi and ii — >■ tX2,3 decays. If both top squarks 
decay into a b and a xf, then 2b + 21 + pr events are 
expected. 

For the light slepton case, results will be shown for 
= 390, 440, 500, 550 and 600 GeV (with masses of 
XiJ X2j X3 at the benchmark values in Table |l|. The heavy 
slepton case considered in the next section has mi_^ — 390 
GeV, with chargino and neutralinos similar to Table |l|. 

We note that there is no stringent constraint on the ti 
mass from ii — > bx^ mode for the benchmark point in 
Table [I with = 500 GeV. The ATLAS search [S] 
in finaTstates with two leptons places limits upto a stop 
mass of ^ 450 GeV in the case where the mass difference 
between ti and x^ is 10 GeV. The ATLAS search [1] in 
1 lepton -I- 4 jets -I- final state places bounds for ii 
mass less than ~ 390 GeV. 

We note, moreover, that these bounds assume B{ti ^• 
&xf ) = 100% in Xi W^*'>x\i whereas for the bench- 
mark scenario considered in Table IIIII in the next Sec- 
tion (" heavy slepton" scenario) with m^^ = 390 GeV, the 
branching fraction to bxi is ~ 60%. Thus, the bench- 
mark point for the " heavy slepton" scenario too is beyond 
the ATLAS bounds. 

The ATLAS top squark searches [4] in single 
lepton-|->4 jet final state assume B{ii — > ix?) — 100%. 
For the top squark masses shown in the paper, the ti 
branching fraction into iX2,3 and ixi are ~40% and ~8% 
respectively. Thus, these benchmark points are not ruled 
out by the relevant ATLAS searches. 

ATLAS top squark searches in the fully hadronic final 
state also assume B{ti txi) ~ 100% and for the 
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benchmark points in this paper, the exclusion Umits are 
weaker than the Umits in the semileptonic final state. 
Thus, the benchmark points are not ruled out by the 
fully hadronic searches either. 

The analysis begins with selecting events with the fol- 
lowing selection cuts: 

(i) At least 2 isolated leptons (e or /i) with > 20 
and 10 GeV in |?7| < 2.5, where the isolation is 
defined as Y.p^rf^"^ < 5 GeV with AR = 0.4; 

(ii) At least 2 jets with > 30 GeV in It^I < 2.5; 

(iii) At least 1 6-tagged jet with px > 30 GeV in |?7| < 
2.5; 

{iv) fk > 150 GeV; 
(v) Ht > 100 GeV. 

At this stage, the dominant SM background is tt 
events. OSSF dileptons arising from the X2 decay are 
kinematically correlated and its dilepton invariant mass 
distribution is expected to have an edge given by 



m; 



edge 



(6) 



The OSSF dilepton mass distribution from tt events 
can be modelled by the dilepton distribution of OSDF 
dilepton events (TH]. The OSSF dilepton mass distri- 
bution from supersymmetric combinatoric background 
{i.e., uncorrelated leptonic pairs) can also be modelled by 
OSDF dilepton mass distribution. This leads to adoping 
subtracting OSDF distributiuon from OSSF distribution. 
The "light slepton" benchmark events would arise in an 
excess in OSSF— OSDF dilepton mass distribution. 

The OSDF dilepton mass distributions for the SUSY 
benchmark point in Table|l]along with SM tt + (0-4) jets 
background is shown in shaded histogram in Fig.[l] while 
its OSSF distribution (blank histogram) is overlayed. A 
clear edge is seen at around 63 GeV for 30 fb^^ luminos- 
ity. 

Figure [2] shows the flavor subtracted distributions at 
AM = m^o - m^o = 53, 63, 70, 77 and 100 GeV for 
ii = 500 GeV, Xi = 113 GeV and m^o - m^o. The 
dilepton mass distribution edge for all these mass differ- 
ences except AM = 100 GeV can be seen clearly. For 
AM = 100 GeV, the signal acceptance is lower and the 
dilepton study does not have sensitivity at 30 fb^^ lumi- 
nosity. 

Figure [3] shows the flavor subtracted distribution for 
tt + (0-4) jets background plus signal events for — 
390,440,500,550 and 600 GeV, with m^o = 113 GeV 



and m 



X3 



X2 



175 GeV. The edge of the dilepton 
distribution for rrif mass upto 550 GeV can be distin- 
guished from the background for 30 fb~^ luminosity. 

An excess for each top squark mass case in Fig. [3] is 
evaluated in terms of significances (5) in Table |ll] We 
show significances in cases where at least one of the jets 
is required to be a 6 -jet. Here S = Ns/\^Ns + Nb, where 




100 
M|| (GeV) 



200 



FIG. 1: The dilepton invariant mass distributions for tt+ (0 — 
4) jets background and the benchmark point in Table |l] are 
displayed for 30 fb~^ luminosity. The unshaded histogram 
shows the M°^^^ distribution, while the shaded histogram 
shows the Mj^^^^ distribution, which is fitted with the dot- 
dashed curve. The solid curve shows the subtracted M^^ 
distribution. 
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FIG. 2: The subtracted dilepton invariant mass distribution 
M^o is varied, for ti = 500 GeV and 



A/,?" as AM = M^o 
xS = 113 GeV for 30 fb" 



luminosity. 
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FIG. 3: The subtracted dilepton invariant mass distribution 
A'/;*!'^ as the i\ mass is varied, all other masses remaining at 
the benchmark value in Table [l] for 30 fb~^ luminosity. 
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A^s and A'b are the number of OSSF dilepton events in 
range of 20 GeV < Mu < 70 GeV for signal (S) and 
background (B), respectively. 

Nb is determined by fitting the entire (SUSY plus tt) 
OSDF dilepton distribution to a polynomial function and 
calculating the number of events in 20 GeV < Mu < 70 
GeV. Ns is the the number of events in excess above Nb- 
We find that the significance of the benchmark scenario 



for mr 



500 GeV is above 3 a for 30 fb" 



TABLE II: [Light slepton case] Signal and background cross 



sections in fb for various ti masses with Xi 



113 GeV and 



-,o ,-0 _ -^rj^ GeV. Significances (5) are given at 30 fb ^ 



ti Mass 


Signal Background 


S 


(GeV) 






(>lb) 


390 


7.08 


46.4 


5.3 


440 


6.0 


45.6 


4.6 


500 (benchmark) 


3.90 


45.1 


3.1 


550 


2.60 


44.9 


2.1 


600 


1.70 


44.8 


1.4 



(i) At least 2 isolated leptons with pT > 20 and 10 
GeV in \r]\ < 2.5; 

(ii) At least 2 jets with pT > 30 GeV in |?7| < 2.5; 

(in) At least one 6-tagged jet with pT > 30 GeV in 
M < 2.5; 

(iv) fk > 190 GeV; 

(w) Ht > 180 GeV; 

(vi) 20 < Mu < 70 GeV. 



In Table IV, we list our signal and background at dif- 
ferent stages of cuts and flavor subtraction. The final 
significance at 30 fb~^ is 0.97 if a & jet is required in the 
event sample. Small value of B{Z — ^ II) causes smaller 
significance in the "heavy slepton" case compared to the 
"light slepton" case. 



TABLE IV: [Heavy slepton case] Cross section (fb) for signal 
and background at diff'erent stages of event selection and fla- 
vor subtractions are shown for the benchmark point in Table 

imi 



III. SCENARIOS WITH HEAVY SLEPTONS 

In this section the case where the sleptons are heavier 
than the lightest top squark in the following relations: 

mr > m^o , m~± > m^o (7) 

The possible ti decay modes are: 

ii ^ txl^tllxl {or tqqxl) (8) 

ii -> bxf~^blDx° {orbqq'xi) (9) 

ii "> txi (10) 

The leptons and quarks from the X2 and Xi^ decays are 
through off-shell Z and W bosons. SUSY masses at a 
benchmark point are shown in Table |III[ 

TABLE III: SUSY masses (in GeV) at "heavy slepton" bench- 
mark point. 



Particle Mass B 

U 390 17% {txD, 14% {txD 

7% (tx?), 62% [hxt) 

Xl 174 7% ml) 

Xl 175 7% ml) 

Xt 164 22% (/z/x?) 

X? 112 



The same final state like in the "light slepton" scenario, 
i.e., 2 jets -I- 2 leptons -I- is still the key in this scenario. 
Therefore the same event selection with re-optimized cuts 
on the ^ and Ht variables. 



Event Selection 


till 


tt + jets 


Ni > 2, Nj >2,Nb> 1, 
fk > 190 GeV, Ht > 180 GeV 


2.1 


84.7 


OSSF dileptons with 
20 < M,?^^^ < 70 GeV 


0.70 


13.2 


OSDF dileptons with 
20 < M,?s°^ < 70 GeV 


0.44 


12.8 


OSSF - OSDF dilepton with 
20 < mOSSf-osdp ^ 


0.26 


0.40 



IV. DARK MATTER RELIC DENSITY 

For the spectrum (light slepton case) in Table |l] Xi 
is a mixture of Bino and Higgsino. In this case, one 
has m^o ^ m-j^o and this proximity of masses increases 
the branching ratio to dileptonic final states. The cor- 
rect relic density may be obtained by choosing the rela- 
tive amounts of Bino and Higgsino in Xi correctly. This 
choice depends on the mass splitting between X2 3 and 
Xi and hence the position of the edge in the dilepton in- 
variant mass distribution. For the benchmark point of 
Table|l], x? is 72% Bino and 28% Higgsino, and the rehc 
density is obtained as ^h? = 0.11, as displayed in the 
first row of Table |Vj The significance of the analysis is 
also given. 

We note that constraints on DM direct detection for 
Bino-Higgsino dark matter depends on the sign of /i, as 
shown for example in |26j . The spectrum considered in 
this paper may be obtained for either sign of /i, and in 
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fact, for fJ, < 0, the DM masses considered in the paper 
are not constrained by current experiments. We note 
that for fJ, > 0, the mixed Bino-Higgsino DM in the mass 
range considered here is ruled out by XENONIOO. 

If the mass sphtting between X2 3 and Xi is increased, 
the rehc density drops due to smaller Higgsino compo- 
nent in The correct relic density may be obtained 
in the light slepton case by coannihilation effects, as dis- 
played in the second row of Table [Vj Due to a low px 
lepton, the significance for this case is smaller. 



For the spectrum (heavy slepton case) in Table III 



Xi is again a mixture of Bino and Higgsino, but due to 
heavier slepton, the annihilation of xi through t-channel 
slepton exchange is not open. Consequently, a larger Hig- 
gsino component in Xi is required to fulfil the relic den- 
sity, and hence the mass splitting between X2,3 ^.nd Xi 
is smaller than the light slepton case. The significance is 
presented in the third row of Table |V) 



TABLE V: Significances (5), relic density (Qh^), and Bino- 
Higgsino composition in three representing SUSY bench- 
maark points for 30 fb~^ luminosity. AM denotes the mass 
splitting between X2 and xi- For sll points, m^o = 113 GeV. 



Masses 
(GeV) 


{B,H} Q.h^ S 
(%) 


Comments 


AM = 64 
m,- = 157 
TUi^ = 500 


{72,28} 0.11 3.1 


Bino-Higgsino DM 
(Light slepton. Table 


AM = 160 
mi = 123 
mr = 500 


{96,4} 0.11 0.44 


Mainly Bino DM 
(Coannihilation) 


AM = 62 
mj = 4000 

mf^ = 390 


{67,33} 0.11 0.97 


Bino-Higgsino DM 
(Heavy slepton. Table HIl 



jets -I- 2 leptons -I- |x with OSSF-OSDF shows larger 
significance when there is a slepton between X2 ^i^d Xi 
which is not unnatural in SUSY models.. 



V. CONCLUSION 

In this paper, we show that top squark searches at the 
LHC8 in dilepton final state could be a powerful probe 
of the composition of x? if 'mix ~ GeV. Motivated 
by Naturalness and the relic density constraint, the spe- 
cific case of a well-tempered neutralino that is a mixture 
of Bino and Higgsino is probed using the decay of top 
squarks. A typical benchmark point has m^o ^113 GeV, 
m^o ~ m^o ^175 GeV, and mi ~ 500 GeV. 

In such scenarios, the lightest top squark decays 
prdominantly into t X2 3 ^^'^ ^ Xi ■ Both X2 and X3 decay 
to II x^ (via an intermdiate slepton state) or Z x^, while 
decays into /i^Xi oi' ^Xi • Therefore, the final state 
has dileptons, jets and missing energy with a clear edge 
in the OSSF— OSDF dilepton mass distribution. This is 
different from the most studied scenarios where x? and 
X2 are mostly Bino and Wino, respectively, and ii deacys 
into txi with branching fraction of 100%. 

In "light slepton" scenario where high yield of dilepton 
events is expected, a discovery sensitivity up to 600 
GeV of TTT-t^ with 30 fb~^ of integrated luminosity at 
the LHC8 is expected. If the Higgsino component in 
x\ is reduced, then one needs coannihilation processes 
to satisfy the relic density. In such a case, the of 
leptons becomes lower and the significance of the study 
is decreased. In "heavy slepton" case, dileptons are 
produced from the Z boson decays. Small branching 
fraction of B{Z II) results in decreasing the discovery 
sensitivity, compared to the "light slepton" case. 



If Xi is totally a Bino, x" a Wino, and the Higgsinos 
are heavy, the relic density tends to be high unless the 
sleptons are light, in which case Xi annihilates through 
t-channel slepton exchange (x? also need to be less than 
100 GeV). However, in this case the branching fraction 
of t to X2 is quite low, and a very large luminosity is 
required for this analysis. 

Table |V] is a summary of the significances for three 
DM scenarios in this study. We see that the signal > 2 
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